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X-ray study of HLX1: intermediate- mass black hole or 
foreground neutron star? 
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ABSTRACT 

We re-assess the XMM-Newton and Swift observations of HLX1, to examine the evi- 
dence for its identification as an intermediate-mass black hole. We show that the X-ray 
spectral and timing properties are equally consistent with an intermediate-mass black 
hole in a high state, or with a foreground neutron star with a luminosity ~ a few 10 32 
erg s _1 ~ 10~ 6 LEdd 7 located at a distance of « 1.5-3 kpc. Contrary to previously 
published results, we find that the X-ray spectral change between the two XMM- 
Newton observations of 2004 and 2008 (going from power-law dominated to thermal 
dominated) is not associated with a change in the X-ray luminosity. The thermal com- 
ponent becomes more dominant (and hotter) during the 2009 outburst seen by Swift, 
but in a way that is consistent with either scenario. 
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1 INTRODUCTION 

Several theoretical arguments have been suggested for the 
formation of black holes (BHs) with masses ~ 10 3 -10 4 M Q , 
straddling the gap between stellar and nuclear supermas- 
sive BHs. However, solid observational evidence of their 
existence remains lacking or disputed. The bright, point- 
like X-ray source 2XMM J011028. 1 -460421 (hen c eforth , 
HLX1 for simplicity), discovered by iFarrell et alj (2009), 
has been proposed as the first unambi guous identification 
of an intermediate-mass BH (see also iGodet et al] [2009; 
IWebb et alj|2010l) . The (variable) X -ray emission from this 
source clearly indicates an accreting compact object rather 
than a star. It appears located insi de or projected in front 



of the bulge/halo of the SO galaxy (Farrcll ct al 
cated at a distance of « 91 Mpc (jAfonso et al 



20091) lo- 
120051 ). If 



HLX1 does belong to that galaxy, it reached X-ray lumi- 
nosities fa 10 42 erg s _1 , implyi ng a BH mass > WOOMp, 
from Eddington- li mit arguments ([Farrell et al.|l2009l ). Its op- 
tical counterpart (|Soria et alj |2010| ) is a point-like source 
with R « 24 mag, which implies an X-ray /optical flux ratio 
~ 500-1000. This is consistent with an X-ray binary, but 
rules out a background AGN, for which we would expect 
flux ratios ~ 0.1-10. The brightness and colour of the op- 
tical counterpart are consistent with either a massive glob- 
ular cluster in ES0243— 49 (which may contain an accret- 
ing intermediate-mass BH) or a foreground M star in the 



E-mail: robcrto.soria@mssl.ucl.ac.uk 



Galactic Halo (|Soria et al.ll2010T ). A residual emission line 
consistent wi th Ha redshifted by t he systemic velocity of 
ES 0243 -49 l|Wiersema et al.ll201ol ) seems to provide deci- 
sive support for the intermediate-mass BH interpretation. 
However, the issue is still hotly debated. 

In this paper, we discuss the constraints to the nature 
of HLX1 provided by the XMM-Newton and Swift obser- 
vations, and re-examine its X-ray luminosity and spectral 
properties. In particular, we want to determine whether X- 
ray flux and spectral information are already sufficient to 
rule out the possibility of a foreground neutron star (NS), 
weakly accreting from a low-mass donor star. If that is the 
case, the X-ray properties of HLX1 could be used in the fu- 
ture as a template to identify other intermediate-mass BHs. 



2 X-RAY OBSERVATIONS 

XMM-Newton's European Photon Imaging Camera (EPIC) 
observed HLX1 on 2004 November 23 (ObsID 0204540201: 
serendipitously and ~ 10' off-axis) and on 2008 November 
28 (ObsID 0560180901: target observation, on axis); see Ta- 
ble 1 for a summary of instrument modes and live times. 
Henceforth, we will refer to those observations as XMM1 and 
XMM2. We downloaded the Observation Data Files from 
the public archive, and used the Science Analysis System 
(SAS) version 9.0.0 (xmmsas_20090615) to process and filter 
the event files and extract spectra. We checked that there 
were no background flares in either observation. For XMM1 
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Date 


Instrument 


Mode 


Live time 


2004 Nov 23 


pn 


prime full (thinl) 


18.0 ks 




MOS 


prime full (thinl) 


21.6 ks 


2008 Nov 28 


pn 


prime small (thinl) 


35.3 ks 




MOS 


prime full (thinl) 


49.9 ks 



Table 1. XMM-Newton observation log. 



MOS, we extracted the source spectra from a circular region 
of radius 45"; for XMM1 pn, we used a 25" x 30" ellipse, 
to reduce overlapping with a chip gap and the row of pix- 
els next to it. For XMM2 MOS and pn, we used a circular 
region of radius 30", because the source is on-axis and has 
a narrower point-spread function. We defined suitable back- 
ground regions to avoid chip gaps. We selected single and 
double events (pattern ^ 4 for the pn and pattern ^ 12 for 
the MOS). After building response and ancillary response 
files with the S AS tasks rmfgen and arfgen, we used XSPEC 
|Arnaudlll996l ) Version 12 for spectral fitting. 

We used the XMM2 data for timing analysis (XMM1 
being too short, and affected by the chip gap problem). The 
science modes used for XMM2 give a time resolution of 5.7 
ms for the pn, and 2.6 s for the MOS detectors. We used 
xmmselect to extract pn and MOS source and background 
lightcurves, selecting single and double events in the 0.2- 
12 keV range. We defined the same start and stop time for 
the pn and MOS lightcurves, so that they could be combined 
into a total EPIC lightcurve. This is possible because pn and 
MOS have a similar exposure times ~ 50 ks, although the 
pn live time is only 71% of the exposure time (small window 
mode). Background subtraction, together with corrections 
for various sorts of detector inefficiencies (vignetting, bad 
pixels, dead time, etc.) was performed with the SAS task 
epiclccorr. We then used standard FTOOLS tasks for timing 
analysis. 

In addition, HLX1 has been the target of over 40 Swift 
X-Ray Telescope (XRT) observations since 2008 October; 
see NASA's Heasarc data archive for a detailed logbook. We 
used the on -line XRT data product generator l|Evans et al.l 
2007, 2009) to extract lightcurves and spectra (including 
background and ancillary response files); we selected grade 
0-12 events. We downloaded the suitable spectral response 
file for single and double events in photon-counting mode 
from the latest Calibration Datab ase (2009 D e cemb er 1); 
it is the same response used by iGodet et al. | <|2009h . We 
grouped the Swift spectra into four bands, according to 
count rates (Section 3.3), and fitted the coadded spectra 
of each band with XSPEC Version 12. 



3 X-RAY SPECTRAL MODELLING 
3.1 Choice of models 

We chose four spectral models based on the combination 
of a power-law with a soft thermal component, suitable 
to the high-accretion-state BH and/or the low-accretion- 
state NS scenarios. The first model is power-law plus 
single-temperature blackbody, which gives the simplest phe- 
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Figure 1. XMM-Newton/EPIC spectra from the 2004 observa- 
tion, simulatancously fitted with the neutron star atmosphere 
model zamp plus power-law. EPIC-pn datapoints and x 2 con- 
tributions are plotted in red; MOS1 data in green; MOS2 data in 
blue. See Table 2 for the best-fitting parameters. 
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Figure 2. As in Figure 1, for the 2008 XMM-Newton/EPIC spec- 
tra. See Table 3 for the best-fitting parameters. 



nomenological estimate of the soft component. The second 
model is power-law plus disk-blackbody, suitable to the BH 
scenario; the thermal disk component has a much broader 
spectral shape than the single-temperature blackbody. The 
models for thermal component in the third and fourth model 
are suitable to fit the emission from a weakly-magnetized 
NS hydrogen atmosphere in hydrostatic and local thermo- 
dynamical equilibrium: they are broader and harder than 
a simple blackbody at the same effective temperature, but 
slightly narrower than a standard disk-blackbody. 

Blackbody (b b in XSPEC), disk-blackbody (diskbb: 



Makishima et alj Il986h and NS atmosphere 



Zavlin. Pavlov fc Shibanovl Il99fj| ) models are well known 
from the standard release of XSPEC and do not require 
additional explanations here. Our fourth model is the zamp 
model, whi ch has been implemented in XSPEC as an addi- 
tive table l|Campana. Mereghetti fc Sidolil l l997T ) using the 
spectra computed in ( Zampieri et al.lll995T ). The difference 
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is that, while the nsa model computes the X-ray spectrum 
of a passively cooling NS, the zamp model was developed 
specifically to reproduce the emission from non-magnetized 
NSs accreting at very low rates (10~ 7 < L/L^dd 10 -3 ), 
for example from the intesterllar medium, a molecular 
cloud, or a very low mass stellar donor (the last case may 
be applicable to HLX1). In fact, the X-ray spectra from nsa 
and zamp turn out to be virtually indistinguishable at the 
signal-to-noise level of our data. The main reason for the 
similarity is that, at such low luminosities, the NS atmo- 
sphere develops smooth temperature and density gradients 
in the inner layers where free-free emission-absorption dom- 
inates; those gradients are very similar in the two cases. The 
free- free opacity is a function of frequency: as a result, the 
emerging higher-frequency photons are emitted in deeper 
(hotter) layers; the observed spectrum is a superposition 
of Planckians at different temperatures, with a broader 
plateau around the peak than a simple blackbody. In 
addition, the zamp model predicts a temperature inversion 
in the most external layers, due to accretion, but since 
this region is already optically thin to X-ray photons, it 
does not appreciably contribute to the observed X-ray 
spectrum. The additional power-law component is observed 
in weakly accreting neutron stars at lumino sities ~ 10 32 
ergs -1 , although its origin is still unclear jjonker et al.l 
120041 '). perhaps associated to faint, residual magnetosphcric 
activity. In the high-state BH scenario, the power-law 
component comes from a hot Comptonizing medium that 
reprocesses part of the thermal disk photons. 

For consistency, we have used the same NS mass M = 
1.4M© and true NS radius R = 12.4 km for both the 
zamp and nsa models. The fitting parameter for the nsa 
model is the local (non-redshifted) effective temperature 
T c ft; the effective temperature inferred by a distant ob- 
server is T cfi -(1 - 2GM/R) 5 . The fitting parameter for 
the zamp model is the total isotropic luminosity at in- 
finity scaled to the Eddington luminosity, L/L^dd- This 
can be easily related to the effective temperature, because 
L = 47r_R 2 o-r e 4 ff (1-2GM/R), assuming that the NS is 
isotropically emitting from the whole surface (a scaling area 
factor can easily be introduced when this is not the case). 
The fitting parameter for the bb model is the colour tem- 
perature Tbb seen by the distant observer (that is, red- 
shifted). Thus, when comparing the best-fitting tempera- 
ture from those three models, we need to remember that 



T bb =7(1- 2GM/R) u b T ci 
hardening factor. 



0.827T o ff, where 7 is the 



3.2 XMM-Newton results 

For both XMM1 and XMM2, we fitted pn and MOS spec- 
tra simultaneously, leaving a free normalization constant 
bewteeen the three instruments. A priori, there might have 
been significant discrepancies due to the pn chip gap cutting 
across the source (in XMM1), or to the fact that pn and 
MOS were in different modes (in XMM2). In fact, the pn 
and MOS events processed and extracted with SAS Version 
9.0.0 turned out to be consistent both in spectral shape and 
normalization (within 3%), for both observations. This is not 
the case when pn and MOS event files are processed with 
earlier versions of the SAS (for example, version 6.6.0 was 
used for the pipeline-processed files in the public archives); 
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Figure 3. Swift/XKT lightcurve from 2008 October to 2010 June, 
with our definition of count-rate bands. See Tables 4-6 for the 
best-fitting parameters to the coadded spectra from bands A, B 
and C. 
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Figure 4. Swift/XRT spectra from band A (red) and band C 
(blue) fitted with the neutron star atmosphere model zamp plus 
power-law. The plot highlights the significant spectral variability 
between the two sets of (coadded) observations. See Figure 3 for 
our definition of Swift/XRT count-rate bands, and Tables 4—6 for 
the best-fitting parameters. 



in particular, for those older versions there is an « 10% dis- 
crepancy between pn and MOS below 0.5 keV, which makes 
it difficult to get strong constraints on t he soft thermal com- 
ponent. This may be the reason why iFarrell et all (2009, 
their Fig. 2) chose to ignore all pn data below 0.5 keV in 
XMM2, and MOS2 data below 0.4 keV in XMM1, without 
providing any explanations. 

The first result of our analysis is that the XMM1 spec- 
trum is dominated by a power- law component, while the 
XMM2 spectrum is dominated by the therma l component: 
this is in ge neral agreement with the results of IFarrell et al.l 
(|2009h and iGodet et al] (|2009h , and confirms their find- 
ing that there was a significant spectr al transition betwee n 
XMM1 and XMM2. However, unlike IFarrell et al.1 l|2009l) . 
we find that the XMM1 spectral fit is improved by the addi- 
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First XMM observation 




wabs*(bb- 


Hpo) 


wabs* (diskbb+po) 


wabs*(zamp+po) 


wabs*(nsa+po) 


Parameter 






Value 








N H ,mt 


4.1I2 9 x 1020 cm ~ 2 


4.3± 3 ' \ x 10 20 cm" 2 3.6 


tf l x 10 20 cm" 2 


3.1± B 2 ° g X 10 20 cm" 2 


kT hh 


n -, 01 +0.028 
u - xo -0.027 


keV 










N hh 












kTdbb 






0.17±°;°4 ke V 








N dhh 






12 9+ 46 - 8 








log (L/L Edd ) 








_6 04+ - 49 






■^zara 








2.01^2 x 10 ~ 5 






kT ca 








[0.059tg : g* keV] 


0.059to| 


$ keV 












1 0+32.3 

1 - 8 -1.4 


x 10- 7 


r 


QQ+ 35 
z - 33 -0.35 


2 Q4+ - 38 
z - y -0.43 


9 R7+0-43 


2. 


S1 +0.51 

-0.28 


JVpo 


6.4± 2 ;° x 10~ 5 


6.1± 2 ; 3 x icr 5 


5.4± 2 ; 4 x 10~ 5 


°- U -1.9 


x 10- 5 



/0.3-10 2 - 7 ±o.3 x 10_13 2 - 7 t° .l x 10- 13 2 -7lo.2 x 10-13 2 - 7 ~ii x 10-13 

/o u .3-io 4 - 3 ±a? x 10-13 4 - 4 ±a8 x 10-13 4 - 1 io.e x 10-13 4 - 3 ±a8 x 10-13 



xl 0.99 (174.1/175) 1.00 (174.7/175) 1.00 (174.8/175) 1.00 (175.3/175) 

Table 2. Best- fitting spectral parameters for the 2004 XMM-Newton/EPIC observation. In addition to the intrinsic column density 
listed here, we included a line-of-sight column density of 2 X 10 20 cm -2 . Errors are 90% confidence level for 1 interesting parameter. The 
fitting parameter for the zamp model is (£/^Edd)i but we have also listed the corresponding value of fcT e ff for an easier comparison with 
the other models. 



tion of a soft thermal component (F-test significance ~ 95% 
for each of the four thermal models used in this work). In 
XMM1, the thermal component carries ~ 1/3 of the emit- 
ted luminosity in the 0.3-10 keV band; in XMM2, it ac- 
counts for almost 80% of the emitted X-ray luminosity. It 
is possible to introduce additional parameters, for example 
using low-metallicity or ionized absorbers, to obtain statis- 
tically equivalent fits without a soft thermal component for 
XMM1. However, the price to pay is that a steeper power- 
law slope is required (photon index T > 3, unusual for an 
accreting compact object). Besides, a combination of ther- 
mal and power-law components is seen in all other obser- 
vations of this source (XMM2 and Swift); so, there are no 
compelling reasons for forcing the first spectrum to be fitted 
by a simple power law. 

When both XMM1 and XMM2 are fitted with a ther- 
mal plus power-law model, we find (Tables 2,3) that the un- 
absorbed flux is the same for both observations, /0.3-10 ~ 
4 X 10" 13 erg cm -2 s" 1 (corresponding to a luminosity 
~ 4 x 10 41 erg s -1 at a distance of 9 1 Mpc), in contrast with 
the findings of iFarrell et all (2009). This is important be- 
cause it means that the spectral transition between XMM1 
and XMM2 was more likely due to a change in the rela- 
tive fraction of thermal/non-thermal photon output at con- 
stant luminosity, rather than to a change in the accretion 
rate. There are subsequent Swift observations when the un- 
absorbed flux does increase, reaching ~ 8 x 10 -13 erg cm -2 
s" 1 (Table 4), but their spectral properties look very differ- 
ent from those seen in XMM1 (Section 3.3). 



The colour temperature of the thermal component is 
« 0.13 keV for both XMM1 and XMM2, when fitted with 
a simple blackbody model, while the effective temperature 
from the zamp and nsa models is ~ 0.06 keV. This is ex- 
pected, because of the hardening effect on the photon spec- 
trum emerging from the NS atmosphere. It implies a hard- 
en ing factor ~ 2.7 , whic h is similar to what was calculated 
bv lZampieri et all l| 19951 . their Table 1) for this range of lu- 
minosities (L « 10 _6 L E dd ~ 2 x 10 32 erg s _1 , in the NS 
scenario). 

The best-fitting intrinsic column density is higher in 
XMM1 (« 4 x 10 20 cm" 2 ) than in XMM2 (< 10 20 cm" 2 ), 
which is at least qu alitatively consistent with the results of 
IFarrell et al ] (|2009l ). We suggest that this may not reflect 
a true physical change in the source. Our X-ray spectral 
models of XMM1 are dominated by the phenomenological 
power-law component, which requires intrinsic absorption 
to avoid divergence at soft energies. We have tried replac- 
ing the power-law model with more complex Comptoniza- 
tion models, where the power-law-like component has an 
intrinsic turnover at low energies, around the temperature 
of the seed photons. We obtain equally acceptable fits (for 
example, = 172.0/175 for a comptt model) with similar 
seed-photon temperatures w 0.1 keV and no intrinsic ab- 
sorption. On the other hand, XMM2 and the Swift spectra 
(Section 3.3) are dominated by the thermal component at 
low energies, and do not require artificial addition of intrin- 
sic absorption. 

The next step is to determine whether the thermal 



© 2010 RAS, MNRAS 000, [[]-?? 



Intermediate-mass BH or foreground NS? 5 



Second XMM observation 




wabs*(bb+po) 


wabs*(diskbb+po) wabs*(zamp+po) 


wabs*(nsa+po) 


Parameter 




Value 





N H ,int 


< 0.5 x 10 20 cm" 2 


< 1.3 x 10 20 cm" 2 


< 0.9 x 10 20 cm -2 


0.9log x 10 20 cm" 2 


kT hh 


0.133lg;°™ keV 








N hh 


3.33+[5; 22 x 10- 6 








kT ihh 




0.186ir 00 1 keV 






N dhh 




20.4l» 9 






loer ( Tj i Tjt^^x i 






-6 02+°-°? 
u -0.01 




-^Vz am 






4.67«; 2 * x 10-5 




kT cB 






[0.060lHo? ke V] 


0.055±°°° 2 keV 


N nsSL 








5.5+^ x IO" 7 


r 


2 62+ 13 
— 0.18 


o riq + 0.23 
^.uo_g 22 


1 89+ 017 
l.o»_ 22 


1 oq+0.18 


N po 


3.4+.°;* x io- 5 




l.5+.°; 3 x 10~ 5 


lAt° -j x 10- 5 


/o.3-10 


3.2+°'} x 10~ 13 


3.4+°'J x 1(T 13 


3-4l° J x 10~ 13 


3-4+o'g x 10~ 13 


JO.3-10 


3.8t° i x 10- 13 


4.0t°' 2 x 10" 13 


3.9±°; 3 x 10" 13 


4.2±°-j x 10- 13 



0.94 (356.5/381) 0.93 (352.6/381) 

Table 3. As in Table 1, for the 2008 XMM-Newton/EPIC observation. 



0.92 (351.7/381) 



0.93 (356.2/381) 



emission is better fitted by a single-temperature or multi- 
temperature component. We find that all four models give 
equivalent fits for XMM1 and XMM2. Taking into account 
all five sets of spectral fits (both XMM-Newton observa- 
tions, and the three grouped datasets from Swift), there is a 
hint that a single-temperature blackbody may give a slightly 
worse fit than broader thermal components (Tables 2-6), 
but longer observations will be necessary to test this sugges- 
tion. In any case, there is no statistical difference between 
disk-blackbody models (most suitable to an intermediate- 
mass BH scenario) and NS atmosphere models. X-ray spec- 
troscopy alone cannot rule out either scenario. 



3.3 Swift results 

Individual Swift/XRT observations do not have enough 
counts to allow two-component spectral fits and to provide 
any constraints on the relative contribution and tempera- 
ture of the thermal component. To get around this problem, 
we examined the Swift/XRT lightcurve and grouped the ob- 
servations into four bands, at very high ("A"), high ("B"), 
intermediate ("C") and low ("D") count rates (Figure 3). We 
then coadded the spectra from all the observations in each 
band. Band D is still too faint for two-component spectral 
fitting, and we will not discuss it here. For bands A, B and 
C we used the same four spectral models applied to XMM1 
and XMM2 (Tables 4-6). The total exposure time of the 
coadded band-A spectrum is 19.0 ks; for band B, 22.3 ks; 
for band C, 67.3 ks. 

We found that a thermal component is required for the 
combined spectrum of every band. In fact, the band-A spec- 



trum (corresponding to the outburst peak in 2009 August) is 
consistent with only a thermal component, without a power- 
law, although the upper limit to the power-law normaliza- 
tion is not very constraining. The fractional power-law con- 
tribution becomes more important for the band-B and band- 
C spectra at lower luminosities, mainly because the thermal 
component declines. Spectral parameters and un-absorbed 
flux of the band-C spectrum are very similar to those of 
XMM2, although at lower signal to noise. 

Putting together the Swift and XMM-Newton spectral 
results, we find two possible trends, which will have to be 
tested by longer observations. First, the relative contribu- 
tion of the power-law component seems to decrease at higher 
luminosities (Figure 5). This can be explained in the frame- 
work of BH accretion, if the Comptonizing region (respon- 
sible for the power-law component) collapses to an optically 
thick disk, at high mass accretion rates. In the framework 
of weakly accreting NSs, our zamp and nsa spectral fits sug- 
gest an un-absorbed 0.3-10 keV luminosity varying between 
« 2-5 x 10 32 erg s -1 (Table 7). In this luminosity range, 
accreting NSs are known to have a power-law and a thermal 
X-ray component, with the relative pow er-law contribution 
decreasing as the luminosity increases (| Jonker et all 12004 
their Fig. 5). Thus, the X-ray spectral evolution of HLX1 
also has some apparent similarities with the behaviour of 
weakly accreting NSs. 

The second possible trend is an increase in the tem- 
perature of the thermal component at higher luminosities 
(Figure 6). This is consistent with optically-thick emission 
from a surface of approximately fixed size, which gets hotter 
perhaps as a result of enhanced accretion rate; it rules out 
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Band-A Swift observations 



wabs*(bb+po) wabs*(diskbb+po) wabs*(zamp+po) wabs*(nsa+po) 



Parameter Value 



N B , int. 


(fixed) 


(fixed) 


(fixed) 


(fixed) 


kT hb 


0.18t°;°3 keV 








N bb 


9.1+2 3 X 10~ 6 








kT Abb 




0'283±°;°2i keV 






N dhh 




8.9±H 






loP" ( T i 1 Jj-CA A 1 






-5 07+°' 1 ? 

J ' u ' —0.13 




zam 






25.2+J' \ x 




kT eS 






[0.104lg;g« keV] 


0.0941°;^ keV 


N nsa , 








1.2+.°;° x 10- 7 


r 


o n + 1.2 
-0.9 


2.0 (fixed) 


2.0 (fixed) 


2.0 (fixed) 


N po 


5.9tg 4 x 10~ 5 


< 4.6 X 10~ s 


< 5.1 X 10~ 5 


< 4.4 X 10~ 5 


/o.3-10 


8.6+2;9 x 10- 13 


7.8+°; 2 x 10" 13 


7.7+°'j x HT 13 


7.7+3 3 x 10" 13 


fun 

J 0.3 — 10 


9-7t 2 ,t x 10 ~ 13 


S.9±o \ x icr 13 


S.7±l[ '% x 10~ 13 


8.8loJ x 10~ 13 



\l 0.76 (14.5/19) 0.74 (15.6/21) 0.67 (14.2/21) 0.69 (14.6/21) 



Table 4. Best-fitting spectral parameters for the coadded Swift/XRT observations at the peak of the 2009 August outburst. See Figure 3 
for our definition of band A. We fixed the intrinsic column density to zero (it converges to zero even when left as a free fitting parameter). 
In addition, we included a line-of-sight column density of 2 X 10 20 cm -2 . Errors are 90% confidence level for 1 interesting parameter. 



Band-B Swift observations 



wabs*(bb+po) wabs*(diskbb+po) wabs*(zamp+po) wabs*(nsa+po) 



Parameter Value 





(fixed) 


(fixed) 


(fixed) 


(fixed) 


kT hb 


°- 14 ±0.02 keV 








N bb 


4.9+J' 3 x 10" 6 








kT dbb 




0.194+™^ keV 






AT d bb 




26.5t 2 «; 2 






log (L/L Edd ) 






-5 O9+0-27 
°- yZ -0.17 




Ni zam 






8.0I3 2 x 10~ 5 




kT cR 






[0.064t°;J 6 keV] 


0.060+.°;°,°? keV 


N nBa , 








5.3jl|;| x 10- 7 


r 


9 c + 0.5 


9 9 + I.O 
z - -1.6 


2.0 (fixed) 


2.0 (fixed) 


JVpo 


5.4+J's X 10~ 5 


2 - 7 -2.6 X 10-5 


2.1±i'2 X 10~ 5 


2.0ti 1 x 10_5 


/0.3-10 


5.0^0 ? x KT 13 


5.1± 2 '° x 10- 13 


5.i±g;g x 10- 13 


5.1±|;° x 10- 13 


/un 

J 0.3 — 10 


5.9+.°;* x 10 ~ 13 


6 -°io.2 x 10-13 


6.0^0 3 x HT 13 


6.lio.4 x 10-13 



1.04 (21.8/21) 1.03 (21.7/21) 0.99 (21.8/22) 1.00 (22.0/22) 



Table 5. As in Table 4, for the 5™/t/XRT observations during the decline from the 2009 August outburst, (as defined in Figure 3). 
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Band-C Swift observations 



wabs*(bb+po) wabs*(diskbb+po) wabs*(zamp+po) wabs*(nsa+po) 



Parameter 



Value 





(fixed) 


(fixed) 


(fixed) 


(fixed) 


kT hb 


0.14±°;°? keV 








N hh 


3.oti.3 x 10-6 








kT dhh 




°- 2 ° 4 ±0.O22 ke V 






N dhh 










loP" ( T i 1 Jj-CA A 1 






°' 1 °— 0.26 




zam 






5.2+ \% x IO" 5 




kT eS 






[0.069±g;~g keV] 


0.066lg;°°* keV 


iV nsa 








2.1+o ° x 10 ~ 7 


r 


o n + 1.2 
-0.9 




i.4±S:i 


1 2 +1 1 


N po 


3.4l^ x icr 5 


i.3±;-° x io- 5 


i.i±; : » x io-6 


0.91q g x 10~ 5 


/o.3-10 


3.2+° l x IO" 13 


3.4+0; 3 x IO" 13 


3.4+° g x IO" 13 


3.5+°;5 x IO" 13 


fun 

J 0.3 — 10 


3.6±°-j x 10- 13 


3.9±0' 3 x IO" 13 


3.9I0 3 x 10~ 13 


4.0t 3 x 10 ~ 13 



0.77 (25.6/33) 0.68 (22.5/33) 0.68 (22.6/33) 0.68 (22.4/33) 

Table 6. As in Table 4, for the 5™/t/XRT observations of luminosity band C (as defined in Figure 3). 



L*p==5E32 



XMM1 



Swift-C 



Switt-B 



XMM2 



Switt-A 



5 10 
Unabsorbed 0.3-10 keV Flux (10~ 13 erg s- 1 cm- 2 ) 

Figure 5. Fractional contribution of the power-law component to 
the un-absorbed 0.3-10 keV flux. The values and error bars have 
been calculated with the zamp model, but very similar results are 
obtained with the other three models. The luminosity markers 
are for a distance of 2 kpc (NS scenario). 



Swift-A 



Swift-B 




XMM2_J_ 




Swift- 




/ 





XMM1 



0.15 



0.2 
T in (keV) 



0.3 



Figure 6. Relation between peak colour temperature of the best- 
fitting diskbb model, and total un-absorbed flux in the diskbb 
component. The dashed line is not a fit to the data: it marks the 
location of the L oc T 4 correlation, expected for optically-thick 
thermal emission from a region of constant area at increasing ac- 
cretion rate (inner disk or neutron star surface). Very similar rela- 
tions between thermal luminosity and temperature are obtained 
with the other three spectral models. 



the alternative possibility of a change in the soft X-ray lumi- 
nosity due to an expanding photosphere. But this scenario 
is equally applicable to emission from the disk around a BH, 
or from a NS surface: in both cases, we expect L ~ T 4 ff when 
the size of the X-ray emitting region is fixed. 



4 SHORT-TERM X-RAY VARIABILITY 

The combined EPIC lightcurve for XMM2 does not appear 
particularly remarkable (Figure 7), but it clearly suggests 
some short-term variability. To quantify such variability, we 
used the background-subtracted pn lightcurve binned to a 
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10 4 



2x10" 3x10 4 4x10" 5x10 4 

Time (s) 




4x1 0- 4 6x10-" 
Frequency (Hz) 



Figure 7. Combined pn + MOS lightcurve in the 0.2-12 keV 
range, binned to a time resolution of 500 s. 



Figure 9. Low-frequency section of the power spectral density for 
the combined EPIC lightcurve. The fitted white noise level has 
been subtracted (normalization parameter = —2 in the FTOOLS 
task powspec). 




-4000 



-2000 

Period (s) - 5504 s 



2000 



Figure 8. A search for periodicities in the combined EPIC 
lightcurve suggests the presence of weak, quasi-periodic modu- 
lations with characteristic periods ss 5300—5600 s, and possibly 
other weak features in the 1000—10, 000 s range. 



time resolution of 1 s. We obtain a Kolmogorov- Smirnoff 
probability of constancy « 8 x 10 -3 , and (53 ±5)% rms frac- 
tional variation in excess of the Poisson level. We applied the 
same analysis to a combined EPIC lightcurve binned to 10- 
s intervals, obtaining a Kolmogorov-Smirnoff probability of 
constancy 1.3 x 10" 2 and (34±2)% rms fractional variation. 

We searched for characteristic periods in the pn and 
combined EPIC lightcurves, by folding the data to a range of 
periods (ef search in FTOOLS). We do not find any strong or 
well defined period; however, there is a weak, quasi-periodic 
modulations with characteristic periods ~ 5300-5600 s (Fig- 
ure 8). The power associated with this modulation or range 
of characteristic modulations is small, barely about the noise 
level, as can be seen from the power spectral density plot 
(Figure 9). At higher frequencies, the power spectral den- 
sity is consistent with white noise. 



5 DISCUSSION AND CONCLUSIONS 

We have re-examined the X-ray spectral evolution of HLX1, 
using two XMM-Newton observations from 2004 and 2008, 
and a series of Swift observations over 2008-2010. In gen- 
eral, a two-component model consisting of power law plus 
soft thermal component provides a good fit to all the spec- 
tra. The first XMM observation is dominated by the power- 
law emission (contributing ~ 2/3 of the flux) but a thermal 
component with a colour temperature ~ 0.13 keV is also 
significantly detected. Conversely, at the peak of the Swift 
outburst, the X-ray spectrum is dominated by a thermal 
component with a colour temperature ~ 0.2 keV. 

Although the relative contributions of thermal compo- 
nent and power-law are significantly different between the 
XMM-Newton observations, the un-absorbed flux is sim- 
ilar, and a factor of two lower than previously claimed 
l|Farrell et al. 1 120 09?). We attribute this discrepancy to our 
processing of the EPIC event files with the latest version of 
the SAS, which may provide more accurate results at ener- 
gies < 0.5 keV. With our re-analysis, pn and MOS now give 
identical spectral parameters and normalizations consistent 
within 3%, even below 0.5 keV. As an aside, we note that 
the spectral difference between the t wo XMM-Newton ob ser- 
vations was visually exaggerated in iFarrell et all $200$ ) by 
their choice of plotting two unfolded spectra based on dif- 
ferent models (simple power-law, and disk-blackbody plus 
power- law). 

For the thermal component, there is no statistical dif- 
ference between disk-blackbody models (most suitable to an 
intermediate-mass BH scenario) and NS atmosphere mod- 
els. X-ray spectroscopy alone cannot rule out either scenario. 
The diskbb model is normalized in terms of apparent inner- 
disk radius (assumed to coincide with the radius of the inner- 
most stable circular orbit, at high accretion rates). Normal- 
iz ation values ~ 10-2 (Tables 2, 3) agree with the findings 
of IFarrell et all |2009f ). and would correspond to a BH mass 
~ 1O 4 M0 at the distance of 91 Mpc, even higher than con- 
servatively estimated from the Eddington limit. The small 
differences (less than a factor of 2) between the fitted radii 
over the five sets of observations may be due to the fact that 
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the inner disk is not always extending precisely to the in- 
nermost stable circular orbit, as well as to other effects such 
as the degree of Comptonization and spectral hardening. If 
HLX1 belongs to ES0243— 49, our spectral fits imply emit- 
ted luminosities ~ 4 x 10 41 erg s _1 in the 0.3-10 keV band, 
during the 2004 and 2008 XMM-Newton observations, rising 
to rss 10 42 erg s _1 at the peak of the August 2009 outburst. 

The zamp and nsa models are normalized in terms of 
distance to the source, assuming uniform emission from the 
surface of a IAMq NS with a radius of 12.4 km. Both mod- 
els suggest a distance « 1.5-3 kpc (Table 7). The small 
discrepancies between the fitted distances over the five sets 
of observations may be due to the fact that the emission 
is never perfectly isotropic from the whole NS surface. If 
HLX1 is a weakly accreting NS, its 0.3-10 keV luminosity 
varies between « (2-5) x 10 32 erg s _1 , that is ~ 10 _6 LBdd- 
This requires accretion rates M ~ 1O _13 M0 yr - 1 . At a dis- 
tance of 2 kpc, the optical counterpart found bv lSoria et al.l 
120101 ) would be a main-sequence M star. We also note that 
any direct X-ray emission from the M dwarf donor would be 
several orders of magnitude below our detection limit: the 
maximum X-ray luminosity of such stars is 10 29 erg s ~ 1 
(IJames et al. 2000;; ISchmit't, Fleming fc Giampapai ri995). 

The increase in the temperature and luminosity of the 
thermal component during the 2009 outburst, and the de- 
crease of the fractional contribution of the power- law compo- 
nent, may be consistent with both a disk-dominated, high- 
state BH and with a quiescent low-mass NS X-ray binary 
(in particular, it is typically seen in NS X-ray binaries with 
X-ray luminosities 10 32 < L x < 10 33 erg s" 1 ). In HLX1, the 
power-law component itself is not constant, as we can see 
by comparing XMM1 with XMM2 and with the very faint 
Band-D Swift observations. In we akly accreting NSs, it was 
suggested, based on the sample of lJonker et al.l l|2004h . that 
the power-law component could represent a constant base- 
line luminosity at « 10 32 erg s _1 for most sources, while the 
thermal component increases with flux. But other quiescent 
NS X-ray binaries contradict this picture, because they de- 
cline to luminositie s (including the po wer-law component) 



< 5 x 10 31 erg s" 1 jHeinke et all 120051) o r even < 3 x 10 



erg s -1 (1H 1905+000~ Tjonker et al.ll2007l ). Evidence of flux 
variability in both the thermal and power-la w components 
is sometimes clearly see n in quiescent NSs |Cackett et al.l 
120051 ; IJonker et all 120051 ); but it is still not clear how the 
variability of the two components may be related, because 
we do not know the physical origin of the faint power-law 
component. 

The X-ray outburst and spectral variability seen in 2009 
may be explained as "canonical" state transitions of BH ac- 
cretion. But they can also be interpreted in terms of spo- 
radic accretion onto the surface of a NS with a low-mass 
donor. Material may accumulate at the magnetospheric ra- 
dius, then be sporadically "flushed" down towards the NS, 
either because of changes in the mass transfer rate from the 
donor star, or changes in the alignment of non-dipole com- 
ponents of the NS field. 

The presence of a variable power-law component (dom- 
inant contribution in XMM1, less important in XMM2 and 
even less in the bright Swift state), with a variable pho- 
ton index « 2-3, does not unambiguosly identify BH or NS 
accretion. The power-law emission may come from inverse- 
Compton scattering of the soft X-ray photons by more 



Observation 


fun 
J 0.3 


-10 


G^zamp 






[W erg 


cm - ^ s~ ^) 


(kpc) 


(kpc) 


XMM1 


4.lj 


-1.0 
-0.6 


2 4 +3.0 
z -^-1.3 


2 4+ 26 


Vl\ /Tl\ /TO 

A MM Z 


3.9^ 


-0.3 
-0.1 


1 a+0-1 
L8 -0.1 


i o+0.2 
L3 -0.1 


Swift-A 


8.7^ 


-1.3 
-0.2 


3 6+ 06 


2 n+0.6 


Swift-B 


6 0+ - 3 
D,u -0.3 


1 6+ ' 5 


1 4+0-5 


Swift-C 


o q+0.3 
o,a -0.3 


n c+0.8 
z -°-0.7 


9 9+O.6 
z - z -0.6 



Table 7. Best-fitting un-absorbed fluxes and distance to HLX1 in 
the NS scenario, from the normalization parameters of the zamp 
and nsa models (Tables 2—6). 



energetic electrons in a hot corona, located either above 
the inner accretion disk or the NS surface. Observation- 
ally, the NS X-ray binary Aql X-l shows X-ray spectral 
variability in quiescence, probably due to variable resid- 
ual accretion: it has a soft thermal component with effec- 
tive temperature varying between w 0.11-0.13 keV, and a 
power-law component with photon index v arying between 
« 1.5 and sa 4, or disappeari ng altogether (|Rutledge et al.l 
120021 : ICampana fc Stellal2003h . Cen X-4 is another quiescent 
NS X-ray binaries well fitted by a soft thermal component 
(kT « 0.16 keV) plus power-law (r ~ 2) when its X-ray lu- 
mino s ity was ~ a few 10 32 erg s _1 l|Campana et al. | |2004 
1 19981 ; lAsai et all Il996l ). Other quiescent NS X-ray bina- 
ries with comparable contributions fro m a soft thermal an d 
a power-law component are listed in IJonker et all <|2004h . 
Quiescent low-mass X-ray binary with X-ray luminosities 
~ 10 32 erg s" 1 , a thermal component and a steep (F > 2) 
power-law component wer e observed in the g lobular clus- 
ter 47 Tuc (W37 and X4: iHeinke et all 120051 ) . Short-term 
X-ray variability w as also seen in some observations of Aql 
X-l (at 3 2% rms; iRutledge et ah 1 12002| ) and Cen X-4 (at 
45% rms; Campana et aTT2004 ). In fact, it was suggested 
l|Heinke et all l2003al lbF that the strength of the power-law 
component and the presence of intrinsic short-term variabil- 
ity in quiescent NS low-mass X-ray binaries are two indi- 
cators of continued low-level accretion. The (speculative) 
detection of weak modulations around ~ 5300-5600 s can 
have many explanations and does not uniquely identify an 
intermediate-mass BH. For example, it is also the orbital 
period of an M4 main-sequence star filling its Roche lobe, 
from the well-known period-dens ity relation in binary sys- 
tems l|Frank. King fc Raine|[2002l ). 

In conclusion, our X-ray spectral and timing analysis 
has provided more accurate constraints on the un-absorbed 
flux, degree of variability, relative thermal/non-thermal con- 
tribution, and temperature of the thermal component from 
HLX1. However, X-ray properties alone are not sufficient to 
rule out either of the competing models previously suggested 
for this source (old NS with a low-mass donor in the Galac- 
tic halo, or intermediate-mass BH in ESO 243— 49). This is 
because the thermal component is equally consistent with 
emission from an accretion disk around a BH or from the NS 
surface, and the presence of an additional power-law com- 
ponent does not unambiguously identify the BH scenario, 
either. The X-ray spectral properties are consistent with an 
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intermediate-mass BH in the high or very high state, but 
also with a quiescent NS with low-level (and variable) resid- 
ual accretion. X-ray variability properties are also consistent 
with both scenarios. We suggest that, regardless of the true 
nature of HLX1, its X-ray properties do not yet provide a 
unique observational signature for the identification of the 
new class of intermediate-mass BHs. Therefore, the iden- 
tification of HLX1 as an intermediate-mass BH must rely 
on its proper ties in other bands, fo r example from optical 
spectroscopy (jWiersema et al.ll2010T ). 
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